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ABSTRACT 

In situ measurements of fluctuation spectra and particle distribution func- 
tions have now been carried out throughout interplanetary space. The link between 
these observations is established by theories of wave-particle interaction. Linear 
instability analysis for the actual non-Maxwell ian particle distribution functions 
and an examination of the velocity dependence of microscopic diffusion coeffi- 
cients form the basis of such an investigation. It is described in more detail for 
the short wavelength, ion acoustic like turbulence which is found by linear in- 
stability analysis to correspond to the observed electrostatic fluctuations. Of 
the transport processes associated with these fluctuations, electron heat conduc- 
tion and electron— ion energy transfer are of particular importance for macroscopic 
solar wind expansion. These effects are studied with the aid o^ an anomalous trans- 
port theory. This theory (Dum, 1978 a, b) is based on the dominance of elastic 
scattering of electrons by fluctuations, similar to (enhanced) electron— ion col- 
lisions. It has a much wider range of applicability than classical transport 
theory, which assumes dominance of Coulomb collisions for elastic and inelastic 
scattering. Nevertheless, a substantial fraction of the heat flux is often carried 
by strongly anisotropic electrons of elevated energy (Strahl) which should be 
handled by a nonlocal transport theory, describing the modification of exospheric 
(collisionless) expansion by residual scattering and the enhancement of the ambi- 
polar electric field by anomalous momentum transfer. 


Solar Wind Models and Breakdown of Classical Transport Theory 

The solar wind like most plasmas we encounter is both turbulent and nearly 
collisionless. Macroscopic studies may point to the importance and direction for 
modifications of classical transport theory but^ due to the interdependence of many 
parameters and the uncertainty of coronal boundary conditions, generally allow no 
unique conclusions on the underlying microscopic physics (Hundhausen, 1972, p,.81). 
The rationale and some principal results from microscopic studies that make direct 
use of the detailed in situ measurements available now for t1le interplanetary so- 
lar wind are outlined in the following. 


The need for a modification of classical transport theory as applied to solar 
wind expansion was recognized at an early stage of theoretical development. (For 
a review see Hundhausen, 1 972.) According to one-fluid theory of radial expansion 
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where p is the mass density, p the pressure, u the speed and the last term in the 
bracket is the gravitational potential, the heat flux q provides the major inter- 
nal source term to which external sources due to hydrodynamic waves propagating 
beyond the corona may be added. The heat flux predicted by classical transport 
theory, however, was much too large to fit observations at 1 AU. Various ad hoc 
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reductions have been introduced to allow for a more efficient conversion of ther- 
mal energy into kinetic energy of expansion. Heat flux is dominated by electrons. 
Two fluid models, however, brought even more problems. Not only were electron 
temperatures too high, but ion temperatures were an order of magnitude too low to 
fit observations at 1 AU. Again, the way out appeared to be an ad hoc increase in 
the electron-ion energy transfer rate. Both modifications have been carried out 
in a number of different ways. The introduction of "fudge factors" in the classic- 
al transport relations may be interpreted in terms of an increase in the effec- 
tive collision frequency, typically by a factor of 10-100 over the frequency of 
Coulomb collisions, which arises presumably by scattering from enhanced fluctua- 
tions. Anomalous transport with even larger increases in the effective collision 
frequency is of course by now a well established fact for many space and labora- 
tory plasmas. It was also recognized very early that solar wind expansion may in 
turn provide the free energy source for enhanced fluctuations. In fact, the vari- 
ous two fluid models themselves imply the onset of microinstabilities already at 
heliocentric distances ranging from 3-11 solar radii, using the skewed electron 
distribution function corresponding to classical heat flux (Forslund, 1970; Singer 
and Roxburgh, 1977). 

Fluid models also imply total breakdown of classical collision dominated 
transport at about the same heliocentric distances even without account of insta- 
bilities, especially for rapidly diverging non-radial flow geometries (Durney and 
Pneuman, 1975). Clearly, heat flux must remain limited as the plasma becomes near- 
ly collisionless, 
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where n is the electron density, Ve=(Te/m)^/^ the electron thermal velocity, 

VeTgi a^(Te^/n) the mean free path for collisions of thermal electrons with ions, 
Lx=(V» In Te)*^ the temperature gradient scale along the magnetic field and Kx a 
dimensionless heat conductivity which depends on the ion charge Z as electron- 
electron collisions with collision time Tee ~ ZTei try to restore a Maxwellian 
distribution. For protons, Z=1 , ^ = 3.16 in the normalization of Braginskii (1967) 
for Tel. The classical (first) relation holds for 1 and Oq = 0(1) cor- 

responds to free streaming. Variations of qe and Aei/Lx between 0.1 and 1 are ty- 
pical for the interplanetary solar wind. The shape of the electron distribution is 
strongly non-Maxwell ian in this case, consisting of a nearly isotropic core and a 
hotter low density halo which may be considered collisionless and may be strongly 
anisotropic (e.g. Feldman et al., 1975, 1979; Ogilvie and Scudder, 1978; Pilipp et 
al., 1981). 


The normalized heat flux q^ depends only on the shape of the electron distri- 
bution, or more precisely the isotropic part fo and the 1=1 anisotropy f]/fo 
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where x = mv /2Tp, fj(x) = (3/2) dcosG f(x,0)cos6. The replacement x — >• x in 

(3) gives (v,j /ve A For the electron rest frame to be used here, (v^ /ve ) =0 thus 
implies fj < 0 at lower energies. 


Classical transport theory assumes the dominance of Coulomb collisions, lead- 
ing to a nearly Maxwellian, fo ~ " (2tt)”3/2 e"^, and weakly anisotropic distri- 
bution, with proportional to A^^/Lx or the parameter (Bx/Z) = (8/3) (^)^/^ 

Agi/LT of Spitzer and Harm (1953). As the collision frequency for large speeds de- 
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creases as these assumptions must break dovm at some speed Vq. This is signal- 
ed by an ever increasing anisotropy (see Fig. 1). 



Figure 1. Heat flux limitation for weakly collisional electrons. Shown are (sche- 
matically): (a) the isotropic part (-.-.-) fo of the distribution f(v) exhibiting 
core and halo, (b) the 1 = 1 component (-) fj(v), which determines the return cur- 
rent of low energy electrons (fi ^O) and the heat flux, due mostly to strongly 
anisotropic high energy electrons. Also shown is fj/fM ( » £m Maxwell- 

ian) predicted by classical transport theory e.g. for A^g/Lx = 0.32 (Aee = Agj 
mean free path for thermal electrons, Lj temperature gradient scale), demonstrat- 
ing total breakdown of this local theory. 

A generous condition for the applicability of the classical transport relation 
in (2) may be found from the requirement that 90% of the predicted heat flux (3) 
should arise from the population v <Vg. This gives Xc = 9 for the Spitzer-Harm 
distribution. Assuming that the validity of the theory can be streched to jf /fM|<l 
gives then the remarkable small values, Ag^/Lx < 10“2, <0.03. Because break- 

down IS rapid (fi/%) = 0(x3), even more generous criteria give similar results, 
see Gray and Kilkenny (1980). For the commonly assumed limit of validity, Ae£/Lx=l 
we obtain xq = 1/2 for ifi/fM| 1, corresponding to a population x < Xg which 
carries a negligible portion of the heat flux. 

In contrast to the predictions of classical transport theory, ifj/fol ■< 3 by 
definition of the 1 = 0,1 components for any f(v). For other physical situations, 
e.g. neutron transport and laser heated plasmas where one has the same problem of 
scale lengths comparable to the mean free path, ad hoc microscopic flux limits 
1^1 /^o I % have long been introduced in order to salvage diffusive transport 
theory. The usual limitation. Of ^ 1 simply comes from the requirement f = fg+fi 
cos > 0 and ignores the fact that Ag£/Lx = 0 (1) also implies significant terms 
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1 ^ 2 in a Legendre expansion of f(v, cos6). The level at which f\/^o saturates 
is not a universal constant but determined by the nearly scatter free motion of 
energetic electrons in the global electric and magnetic fields. Collisionless ex- 
pansion from some exospheric base with magnetic field Bq suggests an approximate 
distribution f(v,S) = f(v), G < 0^, where cosfl^ = (1 -B/Bq)^^^ is the asymptote re- 
sulting from conservation of magnetic moment, inv^2 /2B=const. This model yields 
fj/fo = (3/2) (1+cos^) , independent of energy for x > x^- Observed high energy 
electron distributions are similarly restricted in pitch angle (Strahl) with a 
half width that varies but may become as small as 10^ in the interior of a magnet- 
ic sector. Near isotropy in the vicinity of sector boundaries, however, suggests 
strongly enhanced scattering (Filipp et al., 1981). 

Instabilities and Anomalous Transport 

Breakdown of classical transport theory relatively close to the solar corona 
results in observable strongly non-Maxwe Ilian distribution functions at greater 
distances. Exospheric theory, or more correctly a theory that accounts for the 
rapid increase in Coulomb mean free path with speed, A^£(v)ocv^ (Scudder and 
Olbert, 1979) can in principle account for this non-local transport. However, 
model distribution functions related to these theories or actual observations not 
only allow for essentially the same heat flux instabilities as the weakly non- 
Maxwellian electron distribution for classical heat flux (Schulz and Eviatar, 1972; 
Gary et al., 1975), but open up a wealth of new free energy sources such as ani- 
sotropies and ion beams (see e.g. a recent review by Schwartz, 1980). To be sure, 
instability depends on the assumed range of plasma parameters, but as is well known 
small scale turbulence of varying nature and intensity has now been observed 
throughout the interplanetary solar wind (e.g. Gurnett, 1981). 

Associated anomalous transport in the solar wind as for many other plasmas is 
most frequently discussed by determining plasma parameters corresponding to margi- 
nal stability, or more generally by estimating anomalous relaxation rates for these 
parameters. Such a simplified description is suitable for inclusion in global solar 
wind models (Hollweg, 1978). A simple model for heat flux that includes rates for 
core electron deceleration by the global electric and magnetic fields against ac- 
celeration to the proton frame by Coulomb collisions and the phase speed of waves 
by wave-electron interactions, for example, on comparison with observations indi- 
cates the importance of all three effects in regulating heat flux (Feldman et al., 

1979) . 

If phenomenological models indicate the potential importance of anomalous 
transport then evidently the need for deriving new transport relations from a first 
principle theory is even greater. This difficult task not only involves the calcu- 
lation of particle distribution functions from a kinetic equation with given colli- 
sion terms, as in classical transport theory, but also another kinetic equation for 
the evolution of the wave spectrum. Diffusion coefficients with a velocity depend- 
ence that usually is quite different from that of corresponding terms for Coulomb 
collisions, and wave growth rates for the self-consistent highly non-Maxwell ian 
distribution function are basic new elements in these equations. The important 
questions if and how the detailed microscopic description can be reduced system- 
atically to a fluid type dynamical description that is at least tractable for nume- 
rical codes have been examined previously and were answered in the affirmative for 
the interaction of electrons with ion acoustic like fluctuations (Dum, 1978a, b). 

The key role anomalous transport connected with these fluctuations undoubtely can 
play e.g. in collisionless shocks and laser heated plasmas (e.g. Gray and Kilkenny, 

1980) motivated this research. Because very similar electrostatic fluctuation 
spectra are frequently observed in the solar wind (Gurnett, 1981) their origin and 
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their effect on electron heat flux, electron-ion energy exchange and the global 
electric field should be studied, also with the aim of finding extensions of 
anomalous transport theory. Such a program in no way prejudices the role electro- 
magnetic instabilities can possibly play for these transport processes and certain- 
ly not the role of e.g. anisotropy driven electromagnetic instabilities in shaping 
solar wind ion distributions. 


To see the effect on electrons, ion acoustic like fluctuations may be con- 
sidered an extension of the stable fluctuation spectrum responsible for electron- 
ion collisions with wave number range 1 < ^ = ^e/^o ^ to the range 

Aj)g/ f^< kAjjg < 1 and phase velocities uu/k > (^e I^^bye length, electron gyro- 

radius, bo minimum impact parameter, vi ion thermal velocity). For most electrons 

V > tJu/k, thus scattering is predominantly elastic with diffusion coefficient 0^= 

V v^/2 and frequency (Dum, 1978a) 
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where is the electron plasma frequency, and the average in the second term is 
over the wave spectrum with energy density W = ( 6 E^/8 tt). Wave activity in the 
interplanetary solar wind (Gurnett, 1981) with W/nT^ « 10~7 - 10“5 (peak) complete- 
ly dominates elastic scattering, as A = 1 .24.10^(Tg5 /n)l/2 = \o\\ - 10^^(Tg °K, 
n cm”^). Inelastic scattering with diffusion coefficient D = '%i V 2+Vew<(iu/k)^> 
although comparatively slow, still dominates electron-electron scattering, at least 
for thermal electrons, and tends to flatten the electron distribution. Because of 
the v“3 dependence of , this process is most rapid for low energy electrons, 
taking a time that may be estimated from (Dum, 1978a, Fig. 1) 
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for an initially Maxwellian distribution with temperature 


(5) 


An anomalous transport theory that includes classical transport may be con- 
structed under the assumption that elastic scattering is sufficiently frequent to 
maintain a nearly isotropic electron distribution function against the perturbing 
force^ such as gradients, the electric field etc. (Dum, 1978b). The anisotropic 
part f^of the distribution function becomes a functional of the isotropic part f^, 
with jfVfo i « 1 . In contrast to classical transport theory, £q is generally non- 
Maxwellian and is to be determined from a kinetic equation which includes the much 
slower inelastic scattering and nonlocal propagation effects in time or distance. 

If local relaxation by ion acoustic like fluctuation dominates » then fp takes the 
shape fpocexp -(v/vq)^. The structure of the transport relations depends strongly 
on the shape of fp. For example, in addition to a heat flux proportional to the 
temperature gradient, there is for non-Maxwell ian fo also a heat flux, usually in 
the opposite direction, that is proportional to the density gradient. The effective 
collision frequency, in contrast to the Coulomb collision frequency, is not simply 
a number that depends on a few plasma parameters but is proportional to the fluc- 
tuation level and thus is a dynamic quantity connected with wave growth. Also (4) 
was written for simplicity assuming an isotropic wave spectrum. Actual spectra are 
usually strongly anisotropic and the collision frequency depends thus on the pitch 
angle (Dum, 1978b). 


The anomalous transport theory provides the framework for estimating the effects 
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of the observed short wavelength turbulence in the interplanetary solar wind. 

Taking e.g. W/nTe = 10^7, kAoe = 1/2, n = 17 cm"3 gives Tew = (tt/ 8)1/2 3/\few(ve) = 
25 sec from (4) and Tq ^ 100 Tew from (5). The time Tq for flattening may be still 
longer as the average fluctuation level required by (5) is reduced from the peak 
values by the fact that wave activity occurs in form of many bursts that indivi“ 
dually may be as short as a fraction of a second (Gurnett, 1981). The observed 
shape of fo (Fig. 1) also indicates that the effect of the turbulence on fo is 
nonlocal, except perhaps at very low energies where fo cannot be measured due to 
contamination with photoelectrons and distortions by the spacecraft potential. 
Effective collision times during wave activity are much shorter than the Coulomb 
collision time T^g - Tei (Z=l) = l.SxlO^ sec for Te = 3.6 x1q3 or (parameters cor- 
respond to Fig. 3-5 of Dum et al . , 1980). Coulomb collisions may still play a role 
in regions of velocity space where the resonance condition uHk*y. = 0 is satisfied 

only for a small part of the wave spectrum and obviously during periods of no wave 
activity. 

An estimate for the applicability of the anomalous transport theory may be 
obtained in a similar manner as described above for classical transport. The ano- 
malous theory obviously has a much wider range of applicability because it does 
not demand that fo is Maxwellian and the effective mean free path Aew = Ve/Tew 
for scattering from waves is generally much shorter than for Coulomb collisions 
(A^e - 2 AU for the example given above). Still, the theory assumes dominance of 
elastic scattering and for the v“3 dependence in (4) also must break down at some 
speed, Vc, which, however, is larger than in the case of scattering solely by Cou- 
lomb collisions. 

The detailed measurements of distribution functions available for the solar 
wind are of particular value for an investigation of the physical mechanism of 
heat flux limitation. In order to test the anomalous transport theory, the ratios 
of effective mean free path to the gradient scales of temperature and density as 
well as the ratio of interplanetary electric field to effective runaway field were 
determined from fits of measured and predicted distributions for v "^Vo because 
these values are not known sufficiently well. They come out to be quite reasonable, 
however, also in relation to the prevailing fluctuation levels. For the highest 
measured fluctuation levels among the cases studied so far, using Helios measure- 
ments (Dum et al., to be published), the 1 = 1 component of the distribution func- 
tion predicted by anomalous transport theory for measured fo can be fitted with 
observations in a large velocity range v '*^vc» thus describes almost the entire 
heat flux. For more quiet conditions and large heat fluxes, a substantial fraction 
of the heat flux is carried by energetic electrons v ^ vc> which are strongly ani- 
sotropic. Enhancement of the ambipolar electric field by anomalous momentum trans- 
fer to bulk electrons still affects these electrons, however. 

An instability analysis for waves in the solar wind which directly uses the 
measured distribution functions rather than model distributions has been described 
previously (Dum et al., 1980a, 1981). We now find that most of the observed electro- 
static short wavelength fluctuation may definitely be identified as ion acoustic 
waves. Wave growth and wave characteristics such as uu, i,, depend very strongly on 
the actual detailed distribution functions. This is also demonstrated by the fact 
that we found both weak electron (heat flux) and more recently also ion (beam) 
driven modes, whereas previous calculations with model distribution functions never 
established instability when parameters in these models were not simply assumed but 
determined from fits to actual measurements (Lemons et al . , 1979). The instabili- 
ties we find are very weak, highly variable and limited to relatively narrow ranges 
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)f phase velocity and angles of the wave vector, all consistent with observations 
(Dum et al., 1980b). The short duration of individual wave bursts can be under- 
stood if we close the loop by estimating the time in which e.g. minute changes in 
the electron distribution by isotropization change the sign of the growth rate 
Vk ~ Vke Vki**"^* Lemons et al. (1979). Relaxation oscillation around mar 

gTnal TtabilTty Vk ^ 0 expected if these relaxation effects are combined 

with the perturbing forces that provide the free energy for wave growth. They are 
a common feature of other turbulent collisionless plasmas (Dum, 1981) and imply 
an exchange of energy between electrons and ions, with energy extracted from the 
electrons for Vke > 0» electron (heat flux) driven modes. Rates for energy 

and momentum transfer may be computed from a given wave spectrum once the dielec- 
tric properties are known from the stability analysis with the actual distribution 
functions (e.g. Dum, 1978a, 1981). The rate of momentum transfer (mostly to core 
electrons) is 

^ = Jdk W(k) 2 Im e^(k, u^) k (6) 

where ImCg is the electron contribution to the imaginary part of the dielectric 
constant. In the energy transfer rate Kg, k is replaced by i%. Using the parameters 
given before, Tg/T^ = 7.5, and wave information from Fig. 4 of Dum et al . (1980a) 
gives an acceleration rate du/dt = Rg/nm « 32 km/sec2 as compared to oj/k«64 km/sec 
and u = 392 km/sec. The cooling time is -(3/2)nTg/Kg 3900 sec as compared to 

1 .34 10® sec for energy exchange by Coulomb collisions. It is emphasized again 

that integral effects over many individual wave bursts must be considered, but 
these numbers certainly indicate the potential importance of scattering by ion 
acoustic turbulence. 

In conclusion, we have tried to demonstrate that the recent in situ measure- 
ments of fluctuations and particle distribution functions offer an excellent oppor 
tunity for a test of theories of wave particle interaction and for finding neces- 
sary extensions. 
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